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Context. Nine extrasolar planets with masses between 1 10 and 430Me are known to transit their star. The knowledge of their masses and radii 
allows an estimate of their composition, but uncertainties on equations of state, opacities and possible missing energy sources imply that only 
inaccurate constraints can be derived when considering each planet separately. 

Aims. We seek to better understand the composition of transiting extrasolar planets by considering them as an ensemble, and by comparing the 



Oh 
I 



O " obtained planetary properties to that of the parent stars. 



Methods. We use evolution models and constraints on the stellar ages to derive the mass of heavy elements present in the planets. Possible 



C/2 additional energy sources like tidal dissipation due to an inclined orbit or to downward kinetic energy transport are considered. 

Results. We show that the nine transiting planets discovered so far belong to a quite homogeneous ensemble that is characterized by a mass 
^ of heavy elements that is a relatively steep function of the stellar metallicity, from less than 20 earth masses of heavy elements around solar 
T composition stars, to up to ~ 100 for three times the solar metallicity (the precise values being model-dependant). The correlation is still to 
K> be ascertained however Statistical tests imply a worst-case 1/3 probability of a false positive. 
^ . Conclusions. Together with the observed lack of giant planets in close orbits around metal-poor stars, these results appear to imply that heavy 
■ - - ' elements play a key role in the formation of close-in giant planets. The large masses of heavy elements inferred for planets orbiting metal rich 
stars was not anticipated by planet formation models and shows the need for alternative theories including migration and subsequent collection 
of planetesimals. 
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1. Introduction 

Transiting extrasolar planets are extremely interesting objects 
to study because the coupling of radial velocity and photomet- 
ric measurements allows a determination of both their masses 
and radii, and thus, in principle, a constraint on their composi- 
tion. For that purpose, we need to know as accurately as possi- 
ble the properties of the star-planet systems, and to apply evo- 
lution models using up-to-date input physics. 

Nine transiting extrasolar planets are known to date. In 
Table 1, we summarize their characteristics and that of their 
parent stars on the basis of available literature, including re- 
cently derived metallicities and stellar effective temperatures 
from high-resolution and high signal-to-noise spectra obtained 
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with UVES for the faint star s OGLE-TR-10, 56 111, 113 
and T rES- 1 dSantos et alJ200d) and OGLE-TR- 1 32 jPont et alJ 
'2006^. The ages were derived from isochrone fitting when pos- 
sible. In some cases, lower limits were obtained from lithium 
abundances, the Ca II activit y-age relation, an d an analysis of 
the stellar rotational velocity jMelo et al .'2006'V In the absence 
of useful constraints, we used a conservative upper limit of 
10 Ga because the stars' metallicities and proper motions imply 
that they all belong to the galactic thin disk. 

The radii of close-in extrasolar giant planets (hereafter 
Pegasi planets or Pegasids, after 5 1 Peg b) of a given composi- 
tion will depend on essentially three quantities: their masses, 
ages, and the amounts of flux they receive from their par- 
ent star. The latter are expressed in Table 1 in the form of 
Teq,Q, the effective equilibrium temperature calculated for a 
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Table 1. Characteristics of the transiting Pegasi planets discovered so far 



# 


Name 


Age [Ga] 


[Fe/H] 


7-eff [K] 


Tcq.Q [K] 


Mp[Me] 


Rp [Mm] 


Refs. 


1 


HD209458 


4 


-7 


0.02(3) 


6117(30) 


1487(120) 


210(19) 


97(5) 


[1,2] 


2 


OGLE-TR-56 


2 


-4 


0.25(8) 


6119(60) 


2079(140) 


394(73) 


88(11) 


[3,4,5] 


3 


OGLE-TR-113 


0.7 - 


10 


0.15(10) 


4804(110) 


1340(80) 


429(70) 


77(5) 


[5,6] 


4 


OGLE-TR-132 


0.5 


-2 


0.37(7) 


6210(59) 


1980(100) 


350(38) 


80(!^) 


[6,7,8] 


5 


OGLE-TR-111 


1.1 - 


10 


0.19(7) 


5044(80) 


1033(160) 


168(35) 


71(4) 


[3,5,9] 


6 


OGLE-TR-10 


1.1 


-5 


0.28(10) 


6075(90) 


1578(50) 


200(41) 


102(7) 
81(6) 


[3,4,5,10] 
[11] 


7 


TrES-1 


2 


-6 


0.06(5) 


5226(40) 


1157(140) 


238(22) 


74(1) 


[5,12] 


8 


HD149026 


1.2- 


2.8 


0.36(5) 


6147(50) 


1740(150) 


114(10) 


52(4) 


[13,14] 


9 


HD189733 


0.5- 


10 


-0.03(4) 


5050(50) 


1199(30) 


365(13) 


90(2) 
82.5(23) 


[15] 
[16] 



The numbers in parenthesis represent the uncertainties on the corresponding last digits. 



References: ill .Winn et a 
iBouchv et all j2004); [7] 
j2005h : ri2l lSozzetti et al 
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j2004 : ri3l lSatoetal]j2005h : fUl lCharbonneau et aljiiool) : [ISl lBouchv etai]<2005d) : [16l lBakos et ai]j200d) . 



zero planetary albedo. In terms of these parameters, the sample 
of known transiting Pegasi planets is aheady quite rich. And in- 
deed, evolution calculations of individual planets have yielded 
quite different results. The first one discovered, HD209458b, 
was shown to be anomalou sly large jBodenheimer et al]l200lt 
lOuillot & Showmaiil l2002h . The subsequent six planets ap- 
peared relatively "normal", i.e. fitting the standard evolu- 
tion models within the error bars, wit h the exception of 
OGLE-TR-lOb, another too-large planet tea raffe et aDl2005t 
iLaughl in et al. 2Q05b) - Then, HD 149026b was shown to be 
significantly smaller than expected, requiring the presence 
of a large amount ~ 70Me of heavy elements in its inte- 
rior (Sato et al. 2 0Q1,). The last add ition to that list is that of 
HD189733b IBouchv et alJ l2005bl) . which, like HD209458b 
orbits a star of near-solar metallicity, and is found to be rela- 
tively large. 
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While a radius that is smaller than expected by theoreti- 
cal evolution calculations for a solar-composition giant planet 
can be easily accounted for by the presence of a central dense 
core, or generally more heavy elements in the interior, large 
radii point towards missing physics in the models: an addi- 
tional energy source, or inaccurate physical inputs (equations 
of state, opacities, atmospheric models). It has been proposed 
that the anomalous radius of HD209458b could be explained 
by tidal heating link ed to a small forced eccentricity e ~ 0.03 
jBod enheimer et al] l200ll) . but detailed obser vation indicate 
that the eccentricity is small, e = 0.014 + 0.009 (Laughli n et al] 
l2005ah . and observations of the secondary eclipse imply that 
this would further re quire a chance configuration of the orbit 
jPeming et al.ll20()5l) . Another proposed explanation also in- 
volving tidal dissipation of orbital energy is that the planet may 
be trapped in a Cassin i state with a large orbital inclination 
dWinn & Holm an 2005). Finally, a third possibility that would 
apply to all Pegasi planets is to invoke a downward transport of 
kineti c energy and its dissipation by tides ( Showman & Guillot 
I2OO2 V This last possibility would require the various transiting 
plane ts to have diffe rent core masses to reproduce the observed 
radii (lGuiUol2005l) . 



Fig. 1. Radius anomaly (in units of Mm, or 1000 km) of the 
known transiting Pegasi planets as a function of the metallicity 
of their parent star (see Table 1 for a list of the planets). The ra- 
dius anomaly is defined as the difference between the measured 
radius and the one calculated using a simple evolution model 
for a solar composition gaseous planet. 

The purpose of this article is to test whether these scenarios 
are possible, and to attempt to constrain the masses of heavy 
elements present in these planets. We first compare measured 
radii to the ones obtained by standard evolution models of ex- 
trasolar planets. We derive in Section 3 masses of heavy ele- 
ments and compare them to the metallicities of the parent stars. 
Consequences for our understanding of planet formation are 
then discussed. 

2. Standard models and the radius anomaly 

We first calculate evolution models for all nine planet s based 
on a s imple model assuming a near-solar composition ( iGuillotI 
I2OO5I) . We use a pure hydrogen-helium equation of state 
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dSaumon et alJ 19951) with a helium mass mixing ratio Y - 0.30 
and standard Rosseland opacities jAUardet aljboOll) . The at- 
mospheric boundary condition is calculated for each planet for 
an inc oming stelliu^fl ux averaged over the entire planetary sur- 
face lllroetalJl2005h . We stress that significant uncertainties 
on the atmospheric temperature profile arise from incomplete 
opacity sources, presence or absence of clouds, horizontal in- 
homogeneities, and advective transport of heat. This implies 
that atmospheric temperatures in the literature can vary by up to 
500 K (Barman et al. 2001; Sudarsky et al. 2003; Fortnev et al. 
r2005t llroetalJl2005[ ISeager et alJ l2005l) . For our purposes 
however, this will not matter as long as the temperature of the 
deep atmosphere remains roughly proportional to the equilib- 
rium temperature. For example, in our calculations of atmo- 
spheric profiles, the ratio of the temperature at the 10 bar level 
to the equilibrium temperature is found to be between 1.18 and 
1.28 for all the planets in Table 1. 

We define the radius anomaly as the difference between the 
measured radius and the one obtained by our evolution model 
for the ages inferred in Table 1, and accounting for the fact that 
the m easured radii correspo nd to ~ mbar levels jBurrows et alJ 

E003l .Fortnev et al 200.3: .Iro et alJEoOSi Figure 1 shows the 

radius anomaly thus obtained against the metallicity of the par- 
ent star The error bars on the radius anomaly are equivalent 
Icr error bars calculated on the basis of gaussian distribution 
of radii measurements and on a uniform distribution of the age 
within the limits in Table 1. (Note that for simplicity, in the 
case of OGLE-TR-132b and Tres-lb, we arbitrarily recentered 
the radius measurement to have symetrical error bars). Figure 1 
indicates that the radius anomaly and the metallicity of the par- 
ent star may be anticorrelated. We found no such correlation 
between the radius anomaly and the age, mass, effective tem- 
perature of the star, or with the mass of the planet, and equilib- 
rium temperature. We also found no correlation with the total 
XUV flux received during the planets' lifetime, implying that 
planetary evaporation probably remained limited in magnitude. 

In figure 1 and hereafter, two planets deserve a special 
discussion: OGLE-TR-lOb and HD189733b are characterized 
by two significantly different radius determinations (see ta- 
ble 1). In the case of OGLE-TR-lOb, we adopted the value 
by IBo uchv et al. (2005a), because the value by Holman et al. 
(I2OO5I) remains unconfirmed thus far A special symbol for 
OGLE-TR-lOb is used to remind us of this uncertainty. For 
HD189733b, the most recent value (Bakosetal. 2006) was 
adopted because the photometric analysis appears to be im- 
proved over the old value. In both cases, this corresponds to 
a worst-case situation for the (anti)correlation that we are seek- 
ing. 

The fact that the radius anomaly decreases with increas- 
ing [Fe/H] indicates that the source of the correlation is not 
in the atmosphere because larger amounts of h eavy elements 
would lead to larger atmospheric temperature jFortney et alJ 
I2OO5I) and therefore larger radii. On the other hand, the mag- 
nitude of the anomaly is much too large to be explained by 
assuming that the planets have the same composition as their 
stars. On this basis we now attempt to infer the masses of heavy 
elements present in the planets. 
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Fig. 2. Mass of heavy elements in the planets as a function of 
the metal content of the parent star relative to the Sun. The 
mass of heavy elements required to fit the measured radii is 
calculated on the basis of standard evolution models. Negative 
core masses are required in some cases, implying that some 
significant physical input is missing (see text). Horizontal error 
bars correspond to the Icr errors on the [Fe/H] determination. 
Vertical error bars are a consequence of the uncertainties on the 
measured planetary radii and ages. 

3. Masses of heavy elements 

At this point, the correlation is just indicative however because 
the radius anomaly also depends on the planetary mass. In or- 
der to ascertain its reality and to constrain the planet compo- 
sition/star metallicity relationship, we further calculated evo- 
lution models for planets yyith a ce ntral rock core and a solar- 
composition envelope ('see lGuillotl2005) . The results are shown 
in Figure 2, in which we plot as a function of the metal content 
of the central star, the amount of heavy elements that is required 
to reproduce the observed planetary radii. In most cases, this 
amount is positive, but for HD209458b and OGLE-TR-lOb, the 
large measured radii yield by extrapolation unphysical negative 
values of the mass of heavy elements. 

One possible interpretation is that energy dissipation occurs 
only in HD209458b and OGLE-TR-lOb (e.g. via tidal dissipa- 
tion due to a non-zero inclination of the orbital plane), which 
could explain the negative core masses derived for these plan- 
ets. If we retain this explanation, the remaining 7 planets with 
positive core masses appear to possess amounts of heavy ele- 
ments that are correlated with the stellar metallicity. The corre- 
sponding linear correlation coefficient for these 7 points is then 
relatively high, r = 0.78 (Pearson) or rs - 0.71 (Spearman), 
with a corresponding significance level (no-correlation prob- 
ability) close to 7%. A tough test is however to retry all the 
points within their error bars and calculate an averaged signifi- 
cance level, which we find to be ~ 33%. 

An alternative explanation that we favor is that some 
physics is missing in the treatment of the evolution calcula- 
tion and affects all the planets. One possibiUty could be that 
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Fig. 3. Same as fig. 2, but the mass of heavy elements required 
to fit the measured radii are calculated on the basis of evolution 
models including an additional heat source slowing the cooling 
of the planet. This heat so urce is assumed equal to 0.5 % of the 
incoming stellar heat flux JShowman & Guillol2002l) . 



the hydrogen equation of state severely underestimates the den- 
sity for a given pressure and temperature. It could also be due 
to a strongly underestimated deep opacity that would then ar- 
tificially shorten the planetary cooling and contraction. But 
given that it would require relatively severe modifications to 
these physical inputs, we will only focus here on a third pos- 
sibility: that downward kinetic energy transport and its dissi- 
pation due to tidal efifects in the planetary interior effectively 
leads to an ad ditional energy source that slows the contraction 
of the planets jCiuillot & Showmanl2002l:LShowman & Guillotl 
Quantitatively, we found that adding an additional en- 
ergy source at the planet center equal to 0.5% of the incom- 
ing stellar luminosity was sufficient to solve the problem of 
HD209458b and OGLE-TR-lOb. 

Figure 3 shows the masses of heavy elements that result 
from these evolution models including this additional energy 
source. A pleasing feature of these models is that all plan- 
ets can be explained within the same physical framework, and 
that the planet-star metallicity correlation appears to be rein- 
forced: including all nine planets, the correlation coefficients 
are r = 0.58, rs - 0.50, for a significance level of 17%, which 
is confirmed by permutation and bootstrap tests. Using Monte- 
Carlo retrials of the nine points within their error bars, the mean 
significance level increases to 27%. Removing the questionned 
OGLE-TR-lOb point from the sample, the rank correlation co- 
efficient increases to rs = 0.83, for a false positive probability 
of ~ 1 % and a mean significance level of the Monte-Carlo re- 
trials of 12%. Other calculations (not presented) with a higher 
fraction of dissipated energy ( 1 % of the incoming flux), or with 
an ad-hoc no-helium equation of state also lead qualitatively to 
the same results. 



4. Conclusions 

We have shown that the nine transiting extrasolar planets dis- 
covered so far, including the anomalously large HD209458b 
and the "big core" planet HD 149026b form a relatively ho- 
mogeneous ensemble, whose compositions can be explained 
by the same evolution model for all planets. Alternatively, we 
cannot rule out that the structure of one or several planets was 
affected by particular physical processes (e.g. a non-zero incli- 
nation giving rise to tides, a different orbital evolution. ..etc.). 

We found that the masses of heavy elements in the planets 
appear to be proportional to the metallicities of their parent star. 
This correlation remains to be tested, being still consistent with 
a no-correlation hypothesis at the 1/3 level in the least favorable 
case. If true, it is relatively steep (Mz ~ 60Me(10t^''/'^l - 0.8)) 
and appears to exclude the presence of Pegasi planets around 
stars with a metallicity smaller than that of the Sun by more 
than ~ 0.1 dex, in agreement with the lack of Pegasids (a < 0.1 
AU) around metal-poor stars observed by radial velocimetry. 

The steepness of the relation, and the large values of Mz 
obtained in some cases are surprising and have not been pre- 
dicted by formation models, either in the framework of the 

I ' II 1 1 II 1 

core-instabiHty model ( Ida & Lin 2004; Alibert et al. 2005) or 

by models i nvolving direct gravitational instability of the gas 
llBossl2003l) . The core instability model seems however to pro- 
vide the key ingredients to explain this relation: a faster growth 
of solid planetary embryos with increased metallicity, the sub- 
sequent collection of these during or after the capture of a 
gaseous envelope, and a more likely migration for planets that 
are formed quickly. On the other hand, in the framework of 
gravitational instability models, one needs to assume that the 
gaseous protoplanets are able to collect dust quickly and that 
the efficiency of this process is directly linked to the metallic- 
ity of the disk. However, we would then expect to find planets 
around metal-poor stars which is not observed both by transit 
surveys (although figs. 2 and 3 imply they would be easier to 
detect) and by radial velocimetry. It seems therefore unlikely 
that gravitational instability is responsible for the formation of 
Pegasi planets. 

Another consequence of this work is that the most typical 
Pegasids, 150 M® planets orbiting stars more metal -rich than 
the Sun, will be smaller than expected. Since the transit depth 
scales with the square of the radius, this has tough implica- 
tions for ground-based transit searches. For instance, a 60 Mm 
(0.86 Rj) planet would produce a 0.7% transit in front of a 1 Rq 
star, instead of 1.2% for a 80 Mm (1.15 Rj) planet as typicafly 
assumed, enough to sink it below the detection threshold of all 
present ground-based surveys. And indeed, the planets found 
around the three most metal-rich stars by photometric transit 
searches, OGLE-TR-132, 10, and 56 (see Table 1) are also the 
ones that endure the most irradiation from their parent star: they 
thus have been prevented from contracting below the threshold 
level. This may contribute to explain the relative lack of success 
of photometric transit searches compared to radial velocity sur- 
veys. 

Clearly, although tantalizing, the results presented in this 
study are limited by the relatively small number of known tran- 
siting planets. In the next few years, the discovery of transit- 
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ing extrasolar planets, in particular with the space missions 
COROT (launch 2006) and Kepler (launch 2008), their careful 
characterization coupled to studies of their parent stars should 
greatly improve our understanding of planet formation. 

Note added in proof: The transiting planet recently discov- 
ered by McCullough et al. (2006), XO-lb, is very similar in 
mass, radius and stellar metallicity to HD209458b and further 
strengthens the correlation proposed in this paper. 
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